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Abstract 

The dynamics of the 3 dimensional perfect fluid is equivalent to the motion of 
vortex filaments or "strings". We study the action principle and find that it is 
described by the Hopf term of the nonlinear sigma model. The Poisson bracket 
structure is described by the loop algebra, for example, the Virasoro algebra or the 
analogue of 0(3) current algebra. As a string theory, it is quite different from the 
standard Nambu-Goto string in its coupling to the extrinsic geometry. We also 
analyze briefly the two dimsensional case and give some emphasis on the t«i+oo 
structure. 



1 Introduction 



As Polyakov emphasized in his recent lecture^], the string theory^ may be one of the 
unifying concepts in physics, which is quite essential to the understanding of QCD, quan- 
tum gravity, cosmology or even three dimensional Ising model. The application of string 
theory to these fields are enormously successful, due to its robust mathematical strucure, 
characterized by the Virasoro algebra, the Kac- Moody algebras, BRST symmetry or topo- 
logical methods. In some occations, we could derive exact results just as the consequence 
of their infinite dimensional symmetries. 

Some months ago, Polyakov[@] (See also [^) made a bold proposal that the conformal 
symmetry may determine the critical exponent of the two dimensional turbulence exactly. 
In his argument, only some simple hypothesis together with the conformal invariance gives 
an discrete series of permitted spectrum of the exponent of the turbulence. For three 
dimensional case, MigdalQ proposed to use the loop expansion method of QCD may give 
a nice tool to analyze turbulence. It seems that serious attempt to use string technique 
to the fluid dynamics has just begun. 

In this letter, we would like to propose another "stringy" viewpoints to the Navier- 
Stokes equation. In everyday life, virtually everybody knowns the existence of "string" — 
the vortex filament. It may be observed everywhere, in the motion of tornado, cigaret 
smoke, or in the coffee cup. Usually, this extended object is dissipated after a while. 
However, in the limit where the viscosity vanishes, these filaments can keep their identities 
due to the conservation of circulation. In this case, the field theoretical description of fluid 
reduces to the dynamics of vortex filament. Although the classical motion of this "string" 
has been studied in every detail [Q, we do not seem to have enough knowledge from the 
viewpoint of "string theory" . In this letter, we derive the action principle and the Poisson 
bracket structure for the vortex filament. The kinetic action does not have the standard 
form but is described by the Hopf term of the 0(3) non-linear sigma model[0]. Since 
our action is first order with respect to time, we need to use the canonical formalism for 
the singular system. The algebraic structure, thus obtained, is given by loop algebras, 
such as the Virasoro algebra or an analogure of the 0(3) current algebra. We leave the 
application of our formalism in the future issues. 



2 Two-dimensional vortices and Wi^^o algebra 

As a warm-up, let us think of the dynamics of two-dimensional vortices. It is interesting 
in its own sake since it gives a nice representation of wi+oo algebra. In two dimensions, 
the Navier-Stokes equation takes the following form, 

-wr = "VckOq,, (1) 
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where uj{x) = eapdaUpi^x) is the two dimensional vorticity field and 

Gq,(x) = UJUa{x) — vdaUJ. (2) 

V is the viscocity. 

It is very important that if we start from the delta-function distribution for the vor- 
ticity field, 

N 

i^{x) = ^Tk6{x - Xk{t)), (3) 
fc=i 

it keeps its form under time evolution as long as the viscocity is vanishing. The dynamics 
of the Navier-Stokes equation (field theory equation) thus reduces to the motion of vortex 
"particles", i.e. the equation of motion can be written by Xk{t) alone. To see this, it is 
convenient to combine x and y coordinate hy z = x + iy and z = x — iy. With these 
variables, the equation of motion for the vortex center becomes, 

dzk{t) ^ y dzk{t) r,- 

dt 2Txi .j-f.. zu — dt 2ni zu — 

Although these equations seem quite innocent, it is known that they are far from inte- 
grable. Even for the N=3 case, the motion of the vortices become chaotic. 

These equations reminds us strongly the motion of anyon. Indeed, the Lagrangian of 
the system is in some sense related to the U{1) Chern-Simon term. To derive it, let us 
discuss its Hamiltonian dynamics. The energy functional of the system is given by 



H = ^ f d'^X Ua{x)Ua,{x) 
^ 1 1 



a=l,2 

= --^Y.'^j'^k'^^K^j - - Zk)]. (5) 

The sympletic structure which gives the Hamiltonian equation of motion is 

,, , 1 dU,a) 

l/.*U = -2.i:f-55--). (6) 

The equations of motion Eq.(^) are rewritten as ^ = [zj,H]p^ and ^ = [zj,H]p^. In 
this system, the coordinates z and z are the canonical conjugate with each other. The 
Lagrangian which gives this sympletic structure and the equation of motion is given by, 

i ^ dzi_ 1 

4=1 i^j 

The motion of incompressive fiuid gives the diffeomorphism of two dimensional surface 
which preserve the area element (area preserving diffeomorphism=APD). In the string 
theory, this symmetry is sometimes called wi+oo algebra and plays essential role in several 
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places. Dynamical variables of vortices give a representation of this algebra through 
Poisson bracket. In general, an element of APD is given by the diffeomorphism, 

5z = 5z = (8) 

oz oz 



Since 



-A ^ d^{z,z) _ 2 de{z,z) 

PB~ dz [Zi,(^{Zi,Zi)\p^ - Oij^^ , {'d) 



Zi 1 ^\Z,j , Z-^ 

the generator of the transformation Eq.(^) is given by 



e = - ^ Vje{zi, Zi) = - (fxuj{x)e{x). (10) 

In a sense, it gives the relation between the operator of the quantum mechanics Zi and z 
and the operator of the quantum field theory, uj. This correspondence reminds us of the 
the collective field theory method. The Poisson bracket between uj is 

d d 

[uj{x),uj{y)]pj^ = -2iY,eij — — {uj{y)6{x - y)). (11) 

This is a standard Poisson bracket which defines the APD. 

Recently, Polyakov0 proposed an approach to the two dimensional turbulence based 
on the conformal field theory. In his method, = A~^uj is identified with a primary field 
in a minimal model of CFT. The operator production of the ip field is essential to deal 
with the non-linear term in Navie-Stokes equation. The proposed OPE for ^|J field is that 

Hx[^]~[0] + ---. (12) 

On the right hand side of this OPE, there should not appear the ip field itself. Since ip 
field is the pseudoscalar, it should be true as long as the symmetry of the system is Vir 



X Vir. This statement is somehow in contradiction to our Poisson bracket Eq . (pTD . This 
is possible since our symmetry is APD and the structure constants themselves become 
parity-odd. 

Of course, it is dangerous to think that one may directly apply the Wi+oo symmetry to 
the two dimensional turbulence. In that case, we take the limit when the Reynolds number 
approaches infinity. As a consequence, the viscocity disappears in the final expression. 
However, even if the viscocity is infinitely small, there appears finite amount of outcome. 
It is because the essential nature of the turbulence comes from the balance between the 
dissipation due to the viscosity and the enstrophy flow from outside. As Polyakov|0] 
remarked in his paper, the situation seems to be similar to the anomaly in the quantum 
field theory. 

In this way, the wi+oo symmetry may be broken in the turbulence even if we take the 
limit of vanishing visicocity. It is extremely interesting to investigate how wi+oo is broken 
and possibly the Virasoro structure emerges. 

^ I would like to thank Dr. S. Iso for pointing out this fact. 
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3 Action for vortex filament 



For the three dimensional case, the computation becomes more involved but basical strat- 
egy remains the same. In this case, instead of thinking about point vortices, we need to 
consider the vortex string. The Navier-Stokes equation which is convenient for us is 

^ = -V X G, G^ujxu. (13) 

dt ^ ^ 

The vortex filament is described by the vortex field, 

^i^) = j:r, I da.^^^p^S^'^ {x - Ma,,t)) (14) 

i=l ^ 

a; = V X M 

i=l 



An J dai |f-X- 



One puts these formulae into the Navier-Stokes equation. After some computations, one 
may find the time evolution of the "string field" X{a,t), 

dXj dXj dXj (J^Tk f , dXk Xj — Xk i 

X ^ = ^ X E / d(Tk^ X — ^ . (16) 



daj dt duj y-^ At: J dak \Xj - Xk\^ 
This is equivalent to 

dXj ^ Ffc f dXk Xj — Xk dXi 



dak-=. — X =; h a- 



dt T^TT^ dak \Xj-Xk\^ daj 

r)X 

= u{X,)+a-^. (17) 

a can not be determined from the Navier-Stokes equation alone. This term is present 
because of the freedom of the parametrization of the string aj. 

We claim that the action of the vortex filament is given by, 
S = Sq-H 

So = [dtda,X,.(^x^^^ 



0f 



daj dt 



H ^ - I d^x\u(x) 



= -y^j^k [daJdak^-^] ^ ^ ^ . (18) 
Svr^ ' "J 'J '[da, dak) \X, - Xk\ ^ ' 

As far as we know, this action seems to be new. The second part of the action is identical 
to the Hamiltonian of the system. Together with its variation. 
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and the variation of the first term, Eq.(|T8]) gives equation of motion in the form Eq. ([T6D . 

The first term of the action (Sq) is known as the Hopf term in the theory of 0(3)- 
invariant nonhnear sigma modeL|^ It counts the instanton numbers and describes the 
structure of the 9 vacuum. It is obvious that our 5*0 has also topological nature. Unlike the 
Nambu-Goto string, it does not couple to the two dimensional metric on the world sheet. 
However, it is invariant under the reparametrization of world sheet variables because of 
its Chern-Simon type structure. 

This situation strongly reminds us of the Chern-Simon approach to the braid group [p!o[ 
where Witten has shown that the Jones Polynomial can be described by the correlation 
functions of su{2) Chern-Simon theory. In our case, the fundamental conserved quantities 
of the 3d Navier-Stokes Equation are the energy and the helicity. [|ll|] The latter has definite 
topological meaning since it is expressed by. 



where 



(fxuj{x)u{x) 



2^ 



(20) 



1 



dai I da 



daj 



X 



Xi — Xj 
\Xi-Xj\^ 



(21) 



We remark that is the linking number of two filaments Xi and Xj. Also, the in- 



teresting analogy is the Chern-Simon gravity theory in three dimensions. ||13|| Whereas 
the three dimensional gravity theory is the gauge theory of diffeomorphism of three di- 
mensional space-time, the Poisson bracket of Navier-Stoke equation is volume preserving 



diffeomorphism [H 

Finally, we would like to make a comment on the difference between the vortex motion 
and that of Nambu-Goto string. The issue is their coupling to the extrinsic curvature^. 
To see this, we should first remark that our formulae Eqs.(^^-|T9|) are not well-defined 
due to the divergence in the a integral. For example, the integral which appear in the 
left hand side of Eq.(|r^ diverges logarithmically. We put j = k and write aj = a and 
<Jk = o- + 9. Collecting the most divergent peace gives. 



ax 

'dt 



a=ao 



Att J 2 9 



ax ^ d'^x 

dap dap 
I dX |2 
I dap I 



(22) 



In this form, one recognizes that one may absorb the divergent factor by the redernition 
of time variable, 

t^t = TM, A=llog(L/e), (23) 



This part is a short review. Readers who are famihar with vortex motion can skip the rest of this 
section 
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where L is the length of the vortex and e is the ultraviolet cutoff which may be identified 
with the diameter of vortex filament. In the following, we rewrite t as just t for the 
simplicity. 

The remaining term is directly related to the extrinsic geometry of the curve. Since 

= 1, i.e s gives the 



ax 



the parametrization is arbitrary, we take a special choice of a, 
length of the filament. With this parametrization, we may use the Frenet-Serret formula, 

dX ^ dt , dn - ^ db , ,^ , ^, 

a7"^' d^^^^^^""' = T{s)b - K{s)t, — = -T{x)n. (24) 

One may readily prove that 

ax a^x ( 

^ \ = 4-\- (25) 

|^|2 V / 

This is well-known equation which describes that the vortex filament moves to the bi- 
normal direction with the velocity proportional to its curvature at that point. Very 
interestingly, Hashimoto derived that this equation is exactly solvable. It was proved 
that by changing the variable, 

^{s) = k{s) exp (^i J^' T{s')ds'^ , (26) 

Eq.(p5D becomes the non-linear Schrodinger equation, 

which is integrable. 



4 Canonical formalism 

The action derived in the previous section is the first order in the time variable. Further- 
more, there is the first class constraint because of the reparametrization invariance of cr. 
This is an ideal example to learn the canonical formalism of the singular system. In the 
following, we consider the special case = 1 for the simplicity. The extension to general 
case is straightforward. 

The canonical conjugate of the string field is given by, 

U(a,t)^^ = -e,,,X^^. (28) 



•^The Frenet-Serret formula of the world sheet appeared previously in the context of W-gravity 
theory. The link of vortex-string with W-gravity seems to be a fascinating subject. 
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On the right hand side, there is no time derivative. Hence the definition of momentum 
already gives constraints. Let us denote them as 0p, 

r dX'^ 

<P,{a, t) ^ n,(a) - -e^^^X^-^ ^ 0. (29) 

To see the structure of these constrants, we take the Poisson bracket of them, 

dXP 

[0^(a), M^')]pB = -rw^5(a - a') = M^,{X)6{a - a'). (30) 

Following the general procedure, we add terms proportional to the constraint to the 
Hamiltonian, 

5// = $:/,(x,n)0,. (31) 

p=l 

We need to check whether the time evolution of (f) with respect to new Hamiltonian gives 
rize to any new constraints. The Poisson bracket of the constraint with the extended 
Hamiltonian is, 

3X^ f)X^ 
[0„ H + 5H]pj^ ^ e,^^__«-(x(a)) - Vep^J^^^. (32) 

If we choose = ^m'^(X), the right hand side vanishes weakly. Hence, there are no 
secondary constraints. 

Since the rank of matrix M^j, is two, one linear combination of the constraint fields 
(f)fj_ should be the first class. It corresponds to the reparametrization of the arbitrary 
parameter a. Explicitly, we define, 

dX'' dX^ 
TM.-_^„=-_n„. (33) 

It satisfies the classical version of the Virasoro algebra, 

[T{a),T{a')U = 2T(a')^5(a - a') + ^^S{a - a'). (34) 
Poisson brackets with other fields are given by 

[T{a),X'^{a')U = ^S{a-a'). (35) 
[T{a),U,{a')]^^ = -U,{a)-^6{a-a'). (36) 

[2^(^),0m]pb = -0^(a)^5(a-a'). (37) 

The last equation shows that T{a) is indeed the first class constraint. 

In order to define the Dirac bracket, we need to make gauge fixing of the parametriza- 
tion a. A natural choice is to use the length of the vortex filament as a, 

2 



X 



dX_ 



- 1 f« 0. (38) 
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The Poisson brackets of x with other constraints are 



(39) 



Constraints (f) together with the gauge fixing condition x becomes second class. In general, 
if the second class constraint satisfies Poisson bracket, [ia,ip\pB = Cap, the Dirac 



bracket is given by 



[F, G]^ = [F, GU - E [F^ e.]p^ C-M [e^, G] 



al3 



iPB 



(40) 



a, 13 



In our case, from Eq.(|30D and Eg . (^9]) the matrix G is given by 



GAsicr - a') 



/ 03 -02 bi \ 

—03 Oi 62 

a2 —ai 63 

V -61 -62 -h y 



a, = ~T^6{a - a') 
= -'^£r6{a - a') 



(41) 



The index A stands for (f)^ for y4=l,2,3, and x for A = 4. Because of the constraint x ~ 0? 
the inverse of this matrix becomes very simple. 



/ C3 -C2 di \ 

-C3 ci d2 

C2 -ci (i3 

V -di -d2 -d^ / 



r a. 



^ ^55(- - -') 



(42) 



with 6 is the step function. The Dirac bracket may be computed by using the general 



formula Eq.(40) and, for special case, it gives. 



[X^(a),X'^(a)] 



1 



6ia-a'). 



(43) 



This algebra resembles 0(3) current algebra. However, there appear the derivative of X 
instead of X itself. Furthermore, we should not forget that there is a constraint Eg. (|38D . 



5 Discussion 

There are several issues that we could not discuss well in this letter. Let us make a list 
to illustrate the future issues. 

• Relativistic Formulation: In this letter, our consideration is restricted to the 
non-relativistic situation. Because of this restriction, only one Virasoro algebra could 
appear. If we formulate it relativistically, we expect we have the Virasoro algebras 
both in right and left moving sectors. 
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• Quantization and BRST formalism: Since our algebra Eq.(|38D is slightly dif- 
ferent from the current algebra, the quantization and the BRST analysis seems a 
little bit non-trivial. Anyway, what is "the critical dimension" of our string? 

• Description of Turbulence: As Migdal observed^, there seems to be a relation 
between the viscocity and the plank constant. In this sense, the quantized version 
of our analysis should give some aspects of the turbulence. 

• Relation with the quantum gravity: As we discussed in section 3, there are 
some analogies between three dimensional gravity theory and the Navier-Stokes equa- 
tion. It seems very important to investigate their relation more explicitly. Also, 
in two dimensions, the area preserving diffeomorphism describe some aspect of the 
quantum gravity 0. The appearance of APD in two dimensional Navier Stokes seems 
to establish another link. 



The author would like to thank H. Cateau, T.Eguchi, K. Fujikawa, H. Kawai K. Ogawa 
and Q.-H. Park for their interest on this work. Especially, we are very obliged to S. Iso 
for critical discussions on two-dimensional vortices and Wi+oo algebra and to M. Umeki 
for comments and showing relevant references. 



Note Added: After we submitted this paper, we learned that substantial part of 



this work have been done previously in beautiful papers |]T5[, |jT^. See also [jT^ where 



some material on two dimensional case was discussed. We would like to give our thanks 
to Drs. T. J. Allen, M. Genua, S. Yahikozawa for indicating these references. 
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